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Abstract 

We study the apparent horizon for two boosted black holes in the asymptoti- 
cally de Sitter space-time by solving the initial data on a space with punctures. 
We show that the apparent horizon enclosing both black holes is not formed if 
| the conserved mass of the system (Abbott-Deser mass) is larger than a critical 

C mass. The black hole with too large AD mass therefore cannot be formed in 

the asymptotically de Sitter space-time even though each black hole has any 
inward momentum. We also discuss the dynamical meaning of AD mass by 
q-i examining the electric part of the Weyl tensor (the tidal force) for various 

^ | initial data. 

bX). 
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I. INTRODUCTION 

The inflationary universe scenario is one of the most favorable model to explain the 
present isotropy and homogeneity of the Universe 0. The basic idea of this scenario is that 
the potential energy of a scalar field, which behaves as an effective cosmological constant, 
dominates and causes the de Sitter-like rapid cosmic expansion in the early stage of the 
Universe. Then it seems to be likely that due to this rapid cosmic expansion the initial 
anisotropy and inhomogeneities might be stretched out and the universe becomes homo- 
geneous and isotropic. Such a picture is based on the so-called cosmic no hair conjecture 
(CNHC) which states that "all" space-times with a cosmological constant A approach the 
de Sitter space-time asymptotically [|J. Of course, CNHC is not always true without any 
additional conditions since we know that some inhomogeneities can gravitationally collapse 
into black holes in the space-time with A. Hence the dynamics of the inhomogeneities is 
an important issue to get physical insight into the present homogeneity and isotropy of the 
Universe and there are many researches on this problem ||. 

Recently, we have studied numerically apparent horizons in the initial data in the asymp- 
totically de Sitter space-time [§J - . The results suggest that there is an upper limit on the 
gravitational mass, the Abbott-Deser (hereafter AD) mass of a black hole. The AD mass 
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is the corresponding notion to the ADM mass in the asymptotically flat space-time. Hence 
large inhomogeneities may not collapse into a black hole and, furthermore, large black holes 
may not collide in the asymptotically de Sitter space-time. The same result was obtained by 
the analysis of the Oppenheimer-Snyder model with A ||. Further, dynamical simulations 
for the Brill waves in the asymptotically de Sitter space-time were performed and revealed 
the same results as the spherically symmetric dust collapse ||. It was also shown that under 
some conditions there is an upper limit on the area of the event horizon in the asymptoti- 



cally de Sitter space-time fllQfl . Since, as in the case of the asymptotically flat space-time, 
the area of the black hole does not decrease in the asymptotically de Sitter space-time, the 
black holes with too large total area of the event horizons may not collide and not merge 
if the cosmic censorship hypothesis holds. We note, however, that the relation between the 
AD mass and the event horizon formation is not yet clear in the asymptotically de Sitter 
space-time. 

Although the above analysis show some clear relation between the AD mass and apparent 
horizon formation, there still remain some unanalyzed effects, one of which is initial motion 
of a black hole in two-black hole system. If two black holes have inward velocity, we may 
expect that those black holes might collide each other and will form a single black hole. 
Then, in this paper, we shall make a further study of the AD mass in two-black-hole system. 
Previous work about two-Einstein-Rosen bridge system was focused only on the cases with 
no relative velocity other than the uniform background cosmic expansion ||. Hence, in order 
to see the effect of the relative velocity, we shall investigate the axisymmetric initial data 
of two non-spinning black holes with finite velocity in addition to the background cosmic 



expansion. We solve the initial data on a space with punctures following |[L1||. Then we 



search for the apparent horizon enclosing both black holes in order to get some insight into 
the dynamics of the inhomogeneities in the inflationary universe. We find that such an 
apparent horizon does not appear no matter how close together each black hole is and how 
fast it moves when AD mass of this system is larger than a critical value which agrees with 



that of the Schwarzschild-de Sitter space-time [12]. 

It is worth noting that as pointed out in [13|, AD mass can vanish and furthermore 
can be negative even if inhomogeneities satisfy energy conditions. In this article, we shall 
explicitly show this fact in the case of the initial data containing two black holes. So a 
question arises: "Is gravitational field produced when AD mass vanish?" Then, in order 
to see the dynamical meaning of AD mass, we investigate the asymptotic behavior of the 
electric part of the Weyl tensor which corresponds to the tidal force. In contrast with the 
asymptotically flat space-time, the tidal force comes from the anisotropic velocity as well as 
AD mass. In particular, the contribution from the anisotropic velocity becomes the leading 
term in the asymptotic region when the conserved momentum does not vanish. When AD 
mass vanishes, the tidal force is produced solely by the anisotropic momentum. However, 
even if the dynamical effects of the anisotropic velocity on the large scale inhomogeneities 
might be important, the apparent horizon formation is essentially determined by the AD 
mass and furthermore the critical mass coincides with that of the Schwarzschild-de Sitter 
space-time. This fact suggests that the final state of non-rotating gravitational collapse in 
the asymptotically de Sitter space-time might be the Schwarzschild-de Sitter space-time if 
naked singularities are not formed. Hence, it seems that the inhomogeneities with too large 
AD mass cannot collapse into a single black hole. 
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This paper is organized as follows. In Sec. II, we present initial data of two masses 
with linear momenta in the asymptotically de Sitter space-time and search for the apparent 
horizon, numerically. In Sec. Ill, we investigate the electric part of the Weyl tensor of the 
conformally flat initial data and discuss the effect of AD mass and the linear momentum on 
the electric part of the Weyl tensor. Then Sec. IV is devoted to discussion. We adopt the 



units of c = G — 1. Our notations follow those of MTW [14]. 



II. INITIAL DATA FOR TWO BOOSTED MASSES IN THE ASYMPTOTICALLY 

DE SITTER SPACETIME 

A. Setting up Initial Data 

In order to obtain initial data, we have to solve the Hamiltonian and momentum con- 
straints, 

i3) R-kfki + lK 2 = 6H 2 + 167cp H , (2.1) 
3 

Dj(Ki - 2 -8{K) = 8nJ h (2.2) 

with 

(2-3) 

where is the Ricci scalar of the 3-space. K.{ and K are the traceless part and the trace 
of the extrinsic curvature of the 3-space, respectively. Dj is the covariant derivative with 
respect to the intrinsic metric of 3-space. pu and Ji are, respectively, the energy density 
and the momentum density of the matter fluid measured by the normal line observer. 



In order to solve Eq.(2.1), we follow the York's conformal prescription |15| and assume 
the conformally flat metric, 

df = tp 4 (R, z)f ij dx i dx j , (2.4) 

where /y is the flat Euclidean metric. We set the trace part of the extrinsic curvature K to 

K = -3H. (2.5) 

Note that the condition Eq. ( j2.5| ) turns out to be the Friedmann equation (a /a) 2 = H 2 for 
the scale factor a in the homogeneous and isotropic universe. Hence, the condition (2.5) is 
regarded as the assumption of the uniformly expanding background universe. 
Then Eqs.(2.1) and (2.2) become as 

Aij = -^KjK}ij- 7 - 2np~ H ^- 3 , (2.6) 
Djkl = 8vr Ji, (2.7) 

where = K^ip 2 , p~ H = pu^P 8 , Ji = Ji^P 6 , and L>i is the covariant derivative with respect to 
the flat metric and A is the Laplacian operator of the flat space. The traceless part of the 
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extrinsic curvature K^ is decomposed into the transverse traceless part and the longitudinal 
traceless part 

K id = K T J + (LW)^, (2.8) 

with 

(LW)ij = DiWj + DjWi - p tJ DeW e . (2.9) 

As for the transverse traceless part of the extrinsic curvature, we assume K^ T = 0. Re- 
stricting to the vacuum case, i.e. pn = -h = 0, Eq.(2.7) becomes 

AWi + - Di&Wj = 0. (2.10) 



The analytic monopole solution for W l is given by |16| 

W l = - — {IP 1 + n i P k n k ), (2.11) 
3r 

where n % = x 1 jr. P % is the linear momentum defined by 

P i = j K ij d 2 Sj, (2.12) 

where d 2 Si is the area element of the sphere in the conformal flat space. Just for simplicity, 
we consider only the initial data of two black holes with equal mass and equal magnitude 
of momentum. Since Eq. fl2.10|) is linear, we can easily find the extrinsic curvature for the 
two-black-hole system with linear momenta by superposition as 



K ij = Kg>+Kt\ (2.13) 



where K^ comes from each hole located at (0, 0, ±a). Then K^^ is given as 

#(±)« = JL J>(±> n (±W + p(±)%(±)* - (p - „(±)i n (±)j-)p(±)* n W] , (2.14) 

where n^ 1 and pW 1 are defined by 

n { & i = —{x,y,zTa), (2.15) 
r± 

P (±)i -(0,0,±P), (2.16) 



with r± = yx 2 + y 2 + (z =p a) 2 - Substituting obtained above into Eq.(|2.6|), we get the 
equation for the conformal factor ip. We rewrite the conformal factor following 



^ = -^ + -^ + u = - + u } 2.17 

2r> 2r_ a 
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where Mn is a constant. Then the Hamiltonian constraint becomes 



Am 



-\a 7 kik){l + au) 



(2.18) 



It should be noted that the coefficient in Eq. fl2.18| ) is no longer singular as opposed to 
Eq.( |2.6D because the monopole terms are subtracted by Eq. fl2.17p . The conformal factor 
ip (or u) is obtained by solving iteratively Eq. ( p. 18|) by a finite difference method typically 
with 200 x 200 grid numbers in cylindrical (R, z) coordinate. The boundary conditions are 
given by 



du 
OR 



1 + 



R=0 



M 

2r 
du 

dz 



for r = J x 2 + y 2 + z 2 



+oo, 



0. 



(2.19) 
(2.20) 



Note that 2M agrees with M in the limit P % — > 0. By virtue of the time slicing condition 
Eq.( [2.5|) , the boundary condition Eq.( [2.19|) guarantees that the spacetime is asymptotically 
de Sitter one. In practice, however, numerical solutions cover only a finite region. We 
therefore use the Robin boundary condition at the outer boundary 



du 
dr 



(2.21) 



B. Finding Apparent Horizons 



Now we search for apparent horizons in the initial data prescribed above. A marginally 
trapped surface is a closed 2-surface S where the expansion of future-directed outgoing 
null vectors normal to it vanishes |l7|| . The apparent horizon is defined as the outer 
boundary of a connected component of the trapped region. The apparent horizon is believed 
to agree with the marginally outer trapped surface by the proof in ||17|| . However the proof 
is correct only if the apparent horizon is smooth and hence exactly speaking, the apparent 
horizon might not agree with the marginally outer trapped surface ||18|| . In the situation 
considered here, however, we find a smooth apparent horizon and hence we may regard the 
outer marginally trapped surface as the apparent horizon. 

Let be the outward-pointing spacelike unit normal to S and n M be the unit normal to 



a time slice, then can be written as £ M 



rv 



+ and thus 
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DiS 1 + K ii s i s j 



K = 



on S. For the conformally flat space Eq^^), Eq.( 2.22 ) is rewritten as 

.3 



e = - 



ij 2 (r 2 + r 2 e f/ 2 



+ K ijS l s 3 - K 




+ r e 




— 2r — 4r 



(2.22) 



(2.23) 



where suffices r, 9 denote those partial derivatives. To find an apparent horizon, we solve 
Eq.( [2.23|) as a two-point boundary value problem with boundary conditions rg = at 6 = 
0, 7r/2 following Sasaki et al ||19|| . This method was also used recently by Cook et al. |]20 . 
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Likewise the asymptotically flat space-time, when an apparent horizon forms, there al- 
ways exists an event horizon enclosing it in the asymptotically de Sitter space-time if there 
is no naked singularity and the null convergence condition is satisfied [^TJ. Hence, by inves- 
tigating the existence of the apparent horizon in various initial data, we obtain some insight 
into the evolution of inhomogeneities. 

As in the case of the asymptotically flat space-time, when the separation a between each 
black hole is short enough, the apparent horizon encompassing both black holes appears (see 
Fig.l). The existence of such an apparent horizon means the merging of two black holes and 
we shall search for such an apparent horizon in an asymptotically de Sitter space. 

In order to control the strength of the inhomogeneities, we vary a, and M . However, 
these parameters are not conserved quantities of this system. In the asymptotically de 
Sitter space-time, there are ten conserved quantities associated with ten Killing vectors of 
the background de Sitter space-time. In the initial data considered here, only the AD mass 
associated with the timelike Killing vector of the background de Sitter space-time is the 
non-trivial conserved Killing quantity. (Conserved total linear and angular momentums 
vanish here). Hence, we utilize the AD mass as a measure characterizing the strength of 
inhomogeneities in addition to a and P. 

In general, in conformally flat initial data with the time slicing condition Eq. Q2.5| ), the 
AD mass Mad is expressed as [0,0] 



M AD = M + AM, (2.24) 

with 

M^ — ftfSiFdrf, (2.25) 

AM = — I d 2 Sik l l j . (2.26) 
8n J J 

Here the surface integral is taken over an infinitely large sphere and £* = (x, y, z) is the 
conformal Killing vector orthogonal to the sphere. Although M agrees with the definition 
of the ADM mass and is non-negative, AM can be negative. For example, for the initial 
data that we consider, we obtain from Eq.(2.14) using the Gauss's theorem in the region 
R 3 — {two points} 

H r 

AM = — dvD^KW) + 2HaP = 2HaP, (2.27) 

07T JR 3 -{two points} 

where dv is the volume element in the conformal space and we have made use of Eq. fl2.7| ) 
and the conformal Killing equation (££)«j = 0. Therefore, Mad can be negative for P < 
(inward velocity) if 

M<2Ha\P\. (2.28) 

We shall numerically investigate the critical separation a c such that the apparent horizon 
surrounding black holes disappears for a > a c . In Fig. 2, we display a c jH~ x with respect to 
Mad normalized by M c which is defined by 
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(2.29) 



In fact, M c is the critical mass of the Schwarzschild-de Sitter space-time. The Schwarzschild- 
de Sitter space-time is characterized by Mad-, an d when Mad < M c , it contains a black 
hole. However, in the case of Mad > M c , the Schwarzschild-de Sitter space-time contains 
no black hole. As a test of our numerical code, we searched for an apparent horizon with 
a = and P = 0, which corresponds to Schwarzschild-de Sitter space with Mad = 2M . 
For Mad/M c > 1.007, an apparent horizon is not found. Therefore, our numerical code has 
the accuracy of less than 0.7%. 

In Fig. 2 (a), we set that \P\/M is equal to 1. The white square shows that of P > while 
the black square corresponds to that of P < 0. Fig. 2(b) is the same but for \P\/M = 5. In 
the asymptotically flat case, the critical separation a c with fixed P/Mad is proportional to 
Mad (OK?) and increases monotonically with increasing Mad which coincides with ADM 
mass, since there is no physical scale except for Mad- This behavior means that larger Mad 
produces larger gravity and therefore an apparent horizon is formed even if the separation 
a is large. On the other hand, as shown in Fig. 2, in the case of the asymptotically de Sitter 
space-time, the critical separation a c is not proportional to Mad- There is a maximum 
value near Mad — 0.7 M c and then the critical separation a c decreases with increasing 
Mad- Further, we find that when an apparent horizon forms the inequality Eq. (|2.28 ) is not 
satisfied, that is, Mad is strictly positive. Therefore, we may conclude as far as we have 
examined that if Mad < no apparent horizon encompassing both black holes appears. 

In Fig. 3, (i c /Mad as a function of H/H c is shown, where H c = 1/V27M. We find that 
ci c /Mad drops sharply near H/H c ~ 1. Further, there does not appear the apparent horizon 
enclosing both black holes for all a when Mad > M c (or Had > H c ) within our numerical 
accuracy. This behavior essentially coincides with the two-Einstein-Rosen bridge system 
investigated in ||. When there is a positive cosmological constant, it is impossible for two 
black holes with large Mad to coalesce and to form a single black hole. 



III. THE AD MASS AND LINEAR MOMENTUM AS A SOURCE OF GRAVITY 

In the previous section, we have adopted the AD mass as a measure of inhomogeneities 
and shown its relation to the apparent horizon formation. However, the AD mass can vanish 
and could be negative even if the inhomogeneities due to ordinary matter fields exist. Then, 
we cannot conclude, only from the condition that AD mass vanishes, that such a space-time 
is de Sitter. 

In order to show that Mad can be negative, let us focus on the conformally flat space. 
Then, we will first rewrite AM as 



AM = — 

8tt 



J [CDjKi + itf^LQtf] dv = H J ?Jidv, (3.1) 



where the first equality comes from the Gauss's theorem and in the second one the use 
has been made of the momentum constraint Eq. (|2.2| ) and the conformal Killing equation 
(L£)ij = 0. From Eq.(3.1), we can see that when Jj is negative, AM is also negative. Hence, 
in contrast with the ADM mass in the asymptotically flat space-time, Mad is reduced by 
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the imploding motion of matter fields. Furthermore, in this case, Mad can be negative since 
it is possible to consider arbitrary large H or, in other words, the support of the integrand 
in Eq.(3.1) can be arbitrary large. 

Here we shall estimate the condition to make Mad negative. Assuming the dominant 
energy condition, pjj > (JjJ*) 1 / 2 , we obtain 

\AM\ ~ HL x \ J r \V < HL x p H V, (3.2) 

where L and V are, respectively, the length scale and volume scale of the region in which 
inhomogeneities exist. On the other hand, from Eqs. (|2.6| ) and ( |2.25|) , we obtain 

M = J (p H ^ 3 + -L-KlK)4,- 7 )dv > p H V. (3.3) 

From the above equation, we find the condition of M < \AM\ to be 

L > H-\ (3.4) 

Hence, in order that Mad is negative, the size of the inhomogeneities should be larger than 
the cosmological horizon scale. 

In the case of the asymptotically flat space-time, the ADM mass is the conserved en- 
ergy. Furthermore it can be a source of gravity (tidal force) and therefore it is called the 
gravitational mass. As mentioned, the AD mass is also the conserved quantity and has a 
meaning of energy in the asymptotically de Sitter space-time. However, a question arises 
here: "Does the AD mass produce the gravitational field in the same manner as the ADM 
mass in the asymptotically flat space-time?" In order to answer this question, by analyzing 
the conformally flat initial data, we shall investigate the asymptotic behavior of the electric 
part of the Weyl tensor which corresponds to the tidal force for timelike geodesic normal to 
the spacelike hypersurface considered here. Here we shall consider three cases, i.e., initial 
data in a spherical symmetric space (A), with a spherical source with a linear momentum 
(B), and with two boosted masses of which the total momentum vanishes (C). To obtain 
the "physical" components, we introduce an orthonormal triad of basis vectors e\ a \ as 



4r) = ^q;(x,y,z), (3.5) 



e w = i2 A ! 'j ( xz > y z > ~ x2 ~ y 2 ^ ( 3 - 6 ) 

<* = W* l + v* { ~ y ' x ' 0) - (3 ' 7) 

Using Eq.( [2.5|) , the triad components of the electric part of the Weyl tensor in the vacuum 
region is written as 

E(a)(/3) = C( a) ^) p t^t v = (3) -R( Q )(/3) - 3HK( a )(fl - K\ a) K^) h (3.8) 

where f is the unit vector normal to the hypersurface. 

Since we consider the conformally flat initial data, the asymptotic behavior of the metric 
is given by Eq.( p.l9|) . Hence the triad components of the Ricci tensor ^R^)^ of the 3-space 



asymptotically behave as 



8 



K (r){r) ► IF' ^- y ^ 



( 3 > i2 (e)(e) — ^, (3.10) 



(3) %)M — (3-11) 

and the other components are of order 0(r -4 ). Thus the traceless part of the extrinsic 
curvature K{j gives rise to the difference among the various conformally flat initial data. 

A. Spherically Symmetric Initial Data 

In a spherically symmetric space, Eq. ( [2.7|) is easily integrated as 



S 

K{r){r) = -2K(fi)(fi) = - 2 K(v)(y>) — ►-j, for r ^ +00 (3.12) 



r 

where 

"+00 

„'3 



Sj = 8tt r' 3 J r dr', (3.13) 





and other components of Ki a \m\ vanish. Then, it is easily seen that for r — > +00, Et a )(p\ 
behaves as 

2M 

E( T )(r) = -2E W (o) = -2E {ip){ip) -> — — - HK {r){r) 

= -|(M + |/^) = -^£, (3.14) 



and other components are 0(r 4 ). Eq. ([3.14|) shows that Mad produces the tidal force by 



the same way as ADM mass in the asymptotically flat space-time. 



B. A Spherical Source with Linear Momentum 



The solution of the momentum constraint (|2.7| ) for a single black hole located at the origin 
is obtained by keeping only K^~^ nonzero with a = in Eq. ( 2.13 ). Then P is a conserved 



linear momentum associated with the translational Killing vector of the background de Sitter 
space-time. Then the asymptotic behavior of the extrinsic curvature is given by 

3 

K(r)(r) — ► —Pcos9, (3.15) 
~ 3 

K me) = K( V )( V ) — ► -— Pcos#, (3.16) 

and the other components vanish. We can easily verify that in this case Mad coincides with 
M considering Eqs. (|2. 24 )~( 2. 26]) . However, the asymptotic behavior of £(a)(/3) is different 



from the asymptotically flat case: 
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2 Mad 3H 

£(r)(r)— > ~ 3 —^PcosO, (3.17) 

E m e) = E {v){v) — > -— ^ - |^Pcos0, (3.18) 

and other components vanish. M^d produces the tidal force by the same manner as the 
asymptotically flat space-time. However, in the limit of r — > oo, the leading term of E^m) 
comes from the linear momentum P and depends on the polar angle 6.f\ The relative strength 
between the tidal force due to the momentum P and that due to Mad is given by Hr x 
\P\/Mad- Hence, in order that the tidal force due to the linear momentum dominates, the 
cosmological constant should satisfy \P\/Mad ~ 1 if the inhomogeneity distributes within a 
region smaller than the cosmological horizon scale. On the other hand, if the inhomogeneity 
distributes over the cosmological horizon scale and there is a coherent imploding motion in 
addition to the translational motion, the Mad can vanish. In such a case, the effect of the 
linear momentum on the tidal force also dominates. It should be noted that the tidal force 
is produced by the linear momentum of a matter field even if Mad vanishes. 



C. Initial Data for Two Black Holes with Linear Momenta 

In contrast with the previous case, the total linear momentum vanishes in the case of 

the initial data obtained in Sec. II. From Eqs. fl2.9p and ( 2.11|) , the asymptotic behavior of 
K( a )(p) is 

k {r)[r) — > ^aP(cos29 + 1), (3.19) 

K (em — > -^aP(cos 29 + 3), (3.20) 

K (vKv) ^ -^aP (3 cos 29 + 1), (3.21) 

and the other components are 0(r~ 4 ). Hence Mad is given by 

Mad = M + 2HaP, (3.22) 
and for r — > +oo, is given as 

#(r)(r) — ► - - 1 HaP(3cos20 + 1), (3.23) 



Mad 1 
M AD 1 



E w(o) — " + ^s HaP ( 3 cos2d + 5 )' ( 3 - 24 ) 



E mv) — + —HaP(9cos26 - 1), (3.25) 



1 This effect can be regarded as a kinematical effect due to taking the comoving coordinate system. 
In fact, the linear momentum can be eliminated by a coordinate transformation. We note, however, 
that the dipole term in the tidal force for a test particle cannot be eliminated by the coordinate 
transformation [22]. 
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and the other components are of order 0(r -4 ). Eqs. (|3.23|) - (|3.25|) means that when the 
motion of the inhomogeneities is not isotropic the tidal force depends on their directions 
even if the total momentum vanishes. It is worthy to notice that even if Mad vanishes, the 
tidal force is indeed produced by the anisotropic motion of a matter field as in Sec. III. B. 

IV. DISCUSSION 

We have numerically obtained initial data of two boosted black holes in the asymptot- 
ically de Sitter space-time. We have searched for an apparent horizon enclosing both two 
black holes and have shown the relation between the apparent horizon formation and AD 
mass Mad- Interestingly, when Mad is larger than the critical mass M c of the Schwarzschild- 
de Sitter space-time, the apparent horizon enclosing both black holes does not appear for 
any separation a between each black hole and for any momentum. This result is the same 
as that in the analysis of the Einstein- Rosen bridge system ||. 

Furthermore, in order to understand the dynamical meaning of Mad, we have examined 
the asymptotic behavior of the electric part of the Weyl tensor E^ a )^, which corresponds 
to the tidal force, in three cases of initial data in the asymptotically de Sitter space-time. 
In the spherically symmetric case, Mad produces the tidal force by the same way as the 
asymptotically flat space-time, i.e., Mad/ 7,3 - However, in the case of the initial data of 
a spherical source with linear momentum, the behavior of the tidal force is very different 
from the asymptotically flat case. In this case, although Mad produces a part of the tidal 
force by the same way as the asymptotically flat space-time, the leading term of -E( a )(/3) is 
proportional to r~ 2 , which comes from the linear momentum, and depends on the direction. 
In the case of the initial data obtained in Sec. II, i.e., two boosted black holes, Mad also 
produces the tidal force in the same way as the above examples but the contribution of the 
relative velocity produces the anisotropic dependence of the tidal force on the polar angle 
9. However, in contrast with a spherical source with linear momentum, both contributions 
from Mad and linear momenta on the tidal force have the same asymptotic behavior r -3 
since the total momentum vanishes. 

From these results, it seems that the anisotropic velocity is important for the dynamics 
of inhomogeneities in the scale comparable to the cosmological horizon scale. On the other 
hand, in the analysis of initial data, we have seen that the apparent horizon formation cru- 
cially depends on whether the AD mass is larger than M c or not. This fact does not depend 
on the anisotropic velocity and suggests that the final state of non-rotating gravitational 
collapse in the asymptotically de Sitter space-time is the Schwarzschild-de Sitter space-time 
if the cosmic censorship holds. Inhomogeneities with too large Mad cannot collapse into a 
black hole again if the cosmic censorship holds. 

If there exists a positive cosmological constant, large scale nonspherical gravitational 
collapse may be very different from that of the asymptotically flat space-time due to the 
different behavior of Et a vp\- In order to investigate such an effect on the dynamics of the 
inhomogeneities, we need to solve the Einstein equation numerically and follow the time 
evolution. 
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FIGURE CAPTIONS 



Figure 1: The shapes of apparent horizons near critical separations for P/M = —5 and 
H/2M = 0, 0.1, 0.15 from outside to inside. 

Figure 2: The relation between the critical separation a c / H~ l and AD mass Mad- (a) is 
the case with \P\/M = 1. The white square is for the case with P > and the black square 
corresponds to that of P < 0. (b) is the same as (a) but with \P\/M = 5. The critical 
value looks slightly shifted upward (by 0.7%) due to a numerical error. 

Figure 3: The relation between the critical separation a c /M AD and cosmological constant 
H/H c = JX/A C . 
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